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Abstract 

We investigate the holographic entanglement entropy (HEE) of a strip geometry in four dimen¬ 
sional Q-lattice backgrounds, which exhibit metal-insulator transitions in the dual field theory. 
Remarkably, we find that the HEE always displays a peak in the vicinity of the quantum critical 
points. Our model provides the first direct evidence that the HEE can be used to characterize the 
quantum phase transition (QPT). We also conjecture that the maximization behavior of HEE at 
quantum critical points would be universal in general holographic models. 
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I. INTRODUCTION 


Quantum phase transitions (QPT) are believed to give rise to some of the most interesting 
phenomena in condensed matter physics [1]. Yet a systematic characterization of QPTs 
remains an outstanding open question. Entanglement could provide an important handle 
for understanding QPT. For example, there have been indications that entanglement is 
enhanced at a quantum critical point (QCP) [2-5] (for a review, see [6]). Furthermore, 
change of entanglement patterns may underlie certain QPTs which do not involve symmetry 
breaking or traditional order parameters. 

QPTs are often very difficult to analyze as they naturally occur in strongly correlated 
many-body systems. Entanglement entropy is also notoriously hard to calculate. The 
AdS/CFT correspondence provides powerful tools for both analyzing QPTs and computing 
entanglement entropy [7, 8] in strongly coupled systems. 

In this paper we show that in a class of strongly coupled holographic systems undergoing 
QPTs, entanglement entropy attains a maximum at the corresponding QCPs. Thus entan¬ 
glement entropy can be used as a diagnostic for QPT. Behavior of entanglement entropy 
near various holographic thermal phase transitions were investigated before in [9-12]. 

II. THE HOLOGRAPHIC SETUP AND PHASE DIAGRAM 

Metal-insulator transitions are observed in many condensed matter systems and often 
involve strongly coupled physics. Recently there has been progress in simulating metal- 
insulator transitions in holographic systems in (2-|-l)-dimension [13-16]. Recall that the 
radial direction of an asymptotically Anti de-Sitter (AdS) spacetime can be understood as 
a geometrization of the renormalization group (RG) flow for the dual field theory. The IR 
physics is then encoded in the geometry of the deep interior of the spacetime, to which 
we will refer to as the IR geometry. Thus a metal-insulator transition is realized on the 
gravity side as a geometric transition: as external parameter(s) are dialed, the IR geometry 
deforms from that describing a metallic phase to that corresponding to an insulator. At a 
finite chemical potential, the IR geometry describing a metallic phase is well known, given 
by AdS 2 X [17]. The key to recent development came from construction of the so-called 
Q-lattice models [14] which give rise to new IR geometries dual to insulators. 
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Let us now describe the gravity set up for Q-lattice models. The Lagrangian of the gravity 
dual can be written as 


C = R + 6- - I V<h|2 - 


( 1 ) 


where the AdS length scale is set to unity. $ is a complex held which will be used to introduce 
a lattice structure along one of the spatial directions, say x, by the ansatz $ = e**^*(p 
with if x-independent. A remarkable feature of such an ansatz is that while translational 
symmetry is broken, the gravity equations of motion are still ODEs instead of PDEs, bringing 
great simplihcations. This construction is analogous to the construction of Q-balls which 
is employed to build spherically symmetric solitons in [18], thus is called the holographic 
Q-lattices. 

The full solution to equations of motion of (1) can be written in a form 

\ + Vidx^ + V2d'}/\ , 

V P{^) ) (2) 


A = /i(l — z)adt, $ = e^^^z^ ^0, 

where P{z) = U{1 — z){l + z + z^ — fj?z^/2) and A = 3/2 ± (9/4 + rri^yP. Notice that 
U,Vi,V 2 ,a and 0 are functions of the radial coordinate z only and /i is the chemical po¬ 
tential of the dual held theory. We set the boundary condition for 0 as 0(0) = A, which 
is understood as the lattice amplitude. Then the solutions of the background are specihed 
by three dimensionless parameters, namely the temperature T/p, lattice amplitude A//i^“^, 
and lattice wave number k/fi. For simplicity, we will denote these quantities in short by 
T, A, k through this paper. The metric has an event horizon at = 1 and the spacetime 
boundary is at 2 ; = 0. The Hawking temperature is T = (6 — /i^)17(l)/(87r/i). We will work 
at temperature T = 0.001, which is low enough to observe quantum critical phenomenon. 
In addition, we will set the mass of the scalar held as = —2, with a brief discussion on 
other values of the mass in the end of this paper. 

The phase diagram for the metal-insulator transition over Q-lattice backgrounds (2) can 
be obtained by examining the DC conductivity ctdc of the dual held theory along x direction, 
and is presented in Fig. 1. udc can be computed using standard methods and it can be 
written in terms of various geometric quantities in (2) evaluated at the horizon, i.e. 


o'dc 


V El A;202 ) 

/ 2 = 1 


( 3 ) 
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FIG. 1: The phase diagram over {X,k) plane {X ^ 0,k ^ 0) for the metal-insulator transition in 
Q-lattice background. The upper and the lower region of the plot corresponds to metallic and 
insulating phase, respectively. The HEE with large I is numerically computed through Eq.(6). 
Green triangles represent the peaks of HEE when varying k but fixing A, while red squares are 
peaks of HEE when varying A but fixing k. 

As a practical diagnostic, we classify the dual system into a metallic phase or an insulating 
phase by the temperature dependence of DC conductivity around T = 0.001. More explicitly, 
for the metallic phase we expect StO'dc < 0, while for the insulating phase StO'dc > 0. As 
a result, the QCPs are characterized by c^tO'dc = 0. They form a critical line in k — X plane 
as shown in Fig.l. We stress that this phase structure has little change when decreasing the 
temperature further. 


III. HEE AND NUMERICAL RESULTS 


In the AdS/CFT correspondence the entanglement entropy for a region A is obtained 
from gravity as the area of the minimal surface ja in the bulk geometry which ends at 
dA [7], i.e. 


_ Area(7A) 


(4) 


where is the bulk Newton constant. 
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We now consider the holographic entanglement entropy (HEE) for a strip region in the 
boundary system described by (2). We take strip to stretch in ^/-direction with length 
Ly{-^ cxd) and x-direction with width 21 Ly. While the translational symmetry along x 
direction is broken by <h, the metric components in (2) are still functions of z. It is then 
straightforward to hnd the minimal surface for the strip region, which can be specihed by 
the location of the bottom of the minimal surface in z-direction. We find that satishes 
the equation 


I 



dzz^ 


Vi{z.)V2{z.) 

P{z)Vi{z)W{z,,zy 


( 5 ) 


where W{z^,z) = zlVi{z)V 2 {z) — z'^Vi{zyV 2 {zy. 

From Eq.(4) we hnd that after subtracting the vacuum part the entanglement entropy S 
behaves as 




2/rG7v 



^P{z)Vi{z)W{z,, 


z) 



( 6 ) 


Firstly, we compute the HEE S' as a function of the parameter k and A separately when the 
strip width 21 is hnite but hxed at a relatively large value. Our results are illustrated in Fig. 
2. Interestingly enough, we observe that the HEE displays a pronounced peak in both plots. 
Moreover, we hnd that the location of such turning points is independent of the width of 
the strip when I is relatively large. To show this we plot the shift of the peaks of HEE with 
the width of the strip in Fig. 3. Evidently, the peaks of the HEE converge to a hxed point 
with a dehnite value of k when the width of the strip is becoming large enough. The same 
convergent behavior is also observed in the case when varying the parameter A with k hxed. 
Secondly, we mark the turning points of the HEE with large I in the (A, k) plane of the phase 
diagram as illustrated in Fig.l. Remarkably, we observe that all the turning points of the 
HEE are distributed in the vicinity of the trajectory of the critical points, clearly indicating 
that HEE can be used to characterize the occurrence of QPTs. 

Next, we are interested in the large distance behavior, i.e. I —)■ oo. In this limit it can 
be readily seen from (5) that z* —)■ 1, i.e. the bottom of the minimal surface approaches the 
horizon. It then follows from (6) that S is dominated by the contribution of the part of the 
minimal surface near the horizon. More explicitly, as 1 from (5) and (6), both S and 
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FIG. 2: The left plot is for HEE S' as a function of the lattice wavenumber k with I = 5.56, A = 2, 
while the right plot is for S as a function of the lattice amplitude A with I = 5.44, k = 1.44. 
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EIG. 3: The parameter A is fixed at A = 2 for the above plot. The turning point approaches the 
critical point when the width of the strip is increased from 0.91 to 5.56 uniformly. The maximum 
of each curve is marked as a red dot in figure. 


I are logarithmic divergent, 


S 


VoL 


2J^y 


Vl 


2^Gn V pUB 

I p 


log 


z=l 


+ 


v. 


pUB 


log 


2=1 


I — Z: 


1 — Z: 

+ ■ ■ ■ , 


(7) 


where B = 414 1/2 — — 14^4 p=l + z + z'^ — [i^z^ j^. Thus for large /, 


S — r 14trip + 


with r = 


2Gj^pi? 


( 8 ) 


2=1 


is proportional to the volume 14trip = ILy of the strip, with r the entropy density. Note 
that S is given by the horizon area and thus coincides with the black hole entropy. Thus in 
the large distance limit, S is solely determined by the horizon geometry, which is of course 
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expected. Note that the result is in fact general, applying to generic gravity geometry and 
entangling region of any shape (not just strip). 

In Fig. 4 we give the contour plot of r over the (A, fc)-plane. Again, we see the r achieves 
local maxima precisely at the critical line identihed from the DC conductivity earlier. In 
particular, we show an example in this hgure to locate the maximal value of r along the 
vertical direction (a short green line ) and the horizontal direction (a short red line), which 
are marked by a green dot and a red dot, respectively. In general we hnd the data in 
this method matches rather closely with Fig. 1. In other words, in long distance limit 
entanglement entropy can also be used as a diagnostic for the critical line of QPTs. Moreover, 
in this contour plot we observe that the peaks of HEE form a ridge along the direction of 
the critical line. Thus, when we change a single parameter to locate the position of HEE 
peaks over two or higher dimensional phase diagram, the result obviously depends on the 
direction of observation (or the direction of cutting ridge). Unlike in our current paper, in 
general cases the local HEE peaks obtained in this manner may deviate far from critical 
points in certain directions. Therefore, to avoid this confusion, we may present a more strict 
statement on the relation between HEE and QPTs for a multi-parameter system. That is, 
the peaks of HEE always form a ridge along the critical line, and a pronounced peak can 
be observed at the critical point when computing HEE along the direction perpendicular to 
the critical line over the phase diagram. 

Our result for holographic systems resonates well with earlier observations [2-5] in the 
condensed matter literature regarding entanglement entropy and QCPs. 

Furthermore, we have investigated the behavior of r at lower temperatures. We hnd 
that with the decrease of temperature, the location of the maximum over the phase space 
tends to hx. In metallic phase with hxed parameters the quantity r converges to non-zero 
values. While in insulating phase r tends to be vanishing in zero temperature limit, which 
implies the violation of the volume law of HEE. As a result, the peaks of HEE become 
steep on the insulating side with the decease of the temperature. This phenomenon suggests 
that the maximization behavior of HEE at hnite temperature might be the reflection of the 
discontinuity of HEE at zero temperature. However, as disclosed in [13], in this case the 
quantum phase transition could still be continuous and inhnite order since the free energy 
in the bulk does not exhibit any discontinuous behavior even at zero temperature. Finally, 
it is worthwhile to point out that for insulating phases there could exist a kink (r'(Tc) = 0) 
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FIG. 4: A contour plot of the r, whose values are represented by different colors. Each grey 
trajectory has the same value of r. The blue dashed curve is the critical line where metal-insulator 
transition occurs. 

at Tc ~ 10“® as described in [14]. However, for T < Tc the thermodynamically preferred 
solution still exhibits insulating behavior and r has the singular behavior lini'r-s.o?" —t 0. 
Therefore for even lower temperature, the phase structure illustrated in Figure 1 will be 
stable. 

We conjecture that the discontinuity of HEE at zero temperature may result from the 
artifact of IR geometry AdS 2 - As pointed out in [13, 19], the IR geometry AdS 2 will 
lead to the violation of the third law of thermodynamics due to the non-vanishing entropy 
density at zero temperature, and may therefore be an artifact of large N limit of holography. 
Nevertheless, we intend to argue that the maximization behavior of HEE disclosed in our 
paper may not suffer from this artifact. In addition to the near horizon limit analysis shown 
above. Fig. 3 shows that HEE with hnite I with /T <C 1 also exhibits maximal behavior at 
critical points, which is not dictated by the near horizon geometry since the minimal surface 
can be distant from the horizon. This phenomenon can be understood as that the HEE 
with hnite I actually captures the quantum phase structure at the energy level specihed by 
z^, namely the bottom location of the minimal surface. Furthermore, the peaks of HEE 
converge to hxed values near the QCP when I becomes large (see Fig.l). Therefore, in spite 
of the subtlety of IR geometry AdS 2 , we believe the maximization behavior of HEE at QPTs 
should be genuine and universal in the holographic approach. 
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IV. DISCUSSION 


Throughout this paper we have worked with a scalar mass = —2 because in this 
case we hud the numerical analysis exhibits convergent behavior very well. Nevertheless, 
it is worth to point out that in this setup the AdS 2 BF bound is violated such that new 
phases could be found at lower temperatures [14], Taking this into account, we have also 
performed the calculation for the mass w? = —3/2 which does not violate AdS 2 BF bound. 
We qualitatively obtain the same behavior of the HEE, but with much higher precision and 
more grid points. 

As a summary, in this paper we have computed the HEE of a strip geometry in four 
dimensional Q-lattice background, which exhibits metal-insulator transitions in the dual field 
theory. We have found that the behavior of the HEE in long distance physics is dominantly 
determined by the near horizon geometry of Q-lattice. More importantly, through numerical 
calculation we have demonstrated that the HEE attains a maximum at the corresponding 
QCPs. This fact reveals that HEE can be used as a diagnostic for QPTs indeed. Based 
on what we have observed in Q-lattice setup, we have following remarks and conjectures 
which should be crucial for next investigations on QPT in holographic approach. Firstly, 
we conjecture that the connection between HEE and QPT which is captured by the near 
horizon geometry would be a universal feature for general holographic models. Secondly, 
we further propose a simple but elegant criteria for the occurrence of QPT in holographic 
approach, that would be the existence of a maximum for the area element of black hole 
horizon when changing the parameters of the system in zero temperature limit. Finally, it 
is crucial to explore the scaling behavior of the HEE around QCPs in the zero temperature 
limit, which should be done in a background which contains IR geometry dual to metallic 
phases with vanishing entropy density, rather than AdS 2 with finite entropy density. Our 
investigation on these topics is under progress. 
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